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Abstract A new dynamic control algorithm in order to direct the trajectory of a 
gUder to a pre-assigned target point is proposed. The algorithms runs iteratively and 
the approach to the target point is self-correcting. The algorithm is applicable to any 
non-powered lift-enabled vehicle (glider) travelling in planetary atmospheres. As a 
proof of concept, we have applied the new algorithm to the command and control of 
the trajectory of the Space Shuttle during the Terminal Area Energy Management 
(TAEM) phase. 



1 Introduction 

Space vehicles travel at extreme conditions of speed and acceleration that typically 
do not allow for a "man-in-the loop" approach, forcing, at least partially, automation 
of the flight controls. Thus, automated guidance and control systems are a critical 
component for any re-usable space flight vehicle. 

For example, the implementation of control mechanisms for atmosphere re-entry 
and automatic landing systems used in the Space Shutle focused either on pre- 
programmed manoeuvres following a nominal pre -computed trajectory, or hopping 
across different nominal trajectories whenever the vehicle deviates from an initially 
selected trajectory, [11] and [4]. 

A typical return flight from space has three main phases: 

1) Atmospheric re-entry phase: In this initial re-entry phase the transition from 
spacecraft to aircraft flight mode occurs. The typical altitudes for this phase are 
in the range 120-40 km. 
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2) Glide to the landing site phase, usually referred as Terminal Area Energy Man- 
agement (TAEM), occurring in the altitude range 40-3 km. 

3) Final approach and landing phase, occurring in the altitude range 3-0 km. 

While in the atmosphere re-entry phase, the biggest priority is to ensure that the 
structural constraints of the vehicle are not exceeded; during the TAEM phase, the 
biggest priority is to ensure that the vehicle reaches the Heading Alignment Circle 
(HAC) where preparation for landing is initiated. 

On a typical mission, the TAEM phase begins at the altitude of 25,000 — 
40,000 m at a speed around 2 - 6 M (Mach), and finishes at the HAC at the al- 
titude of 1 , 500 - 3 , 000 m, with a speed of the order of 0.20 M. 

In this paper, we propose a new dynamic control algorithm in order to redirect the 
trajectory of gliders to a pre-assigned target point. This algorithm runs iteratively en- 
abling a self-correcting approach to the HAC and is applicable to any non-powered 
lift-enabled vehicle (glider) travelling in planetary atmospheres. 

This paper is organised as follows. In section 2, we present the equations of mo- 
tion of a ghder and we discuss the approximations we use to define the controllabil- 
ity conditions. In section 3, we briefly discuss the phenomenology of aircraft gliding 
motion, instrumental for the design of a dynamic control strategy. In section 4, we 
derive the dynamic control algorithm, and in section 5 we present realistic simu- 
lations for the Space Shuttle TAEM guidance and control. Finally, in section 6 we 
discuss the main conclusions of the paper. 



2 Gliding motion 

We consider that aircraft gliding motion in a planetary atmosphere is well described 
by a point mass vehicle model under the influence of a gravity field, [3], [5], [9] and 
[2]. In this case, the equations of motion of a gliding aircraft (no thrust forces) are, 

mV — —mg{z)sinY — D{a,Ma) ( x = Vcos^cosy 

mVY ~ —mg{z) CO?, J -\-L{a, Ma) co^H ,< j = V sin cosy (1) 
my;t;cos7 = Lsin;U I z = ysin7 



where m is the aircraft mass, V = ^Jv} + Vy -|- V^-^ is the aircraft speed, y is the flight 
path angle as defined in figures 1 and 2, ji is the bank angle as defined in figure 2c), 
D{a,Ma) and L{a,Ma) are the drag and Uft forces induced by the atmosphere, a 
is the angle of attack and Ma is the Mach number. In general, the Mach number Ma 
is a function of V and z. The function g{z) = goi^E / {Re + z))^ is the gravity accel- 
eration, go = 9.80665 m/s^ is the Earth standard gravitational acceleration constant 
and /?£ = 6.371 x 10^ m is the Earth (or planetary) mean radius. 

In the local reference frame of the aircraft, figure 1,V E (0,°°), 7€ [—n/2, n/2] 
and X G [0,27r]. The bank angle is defined in the interval [—n/2, n/2]. In this 
reference frame, positive values of n correspond to left tums and negative values of 
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H correspond to right turns. As usual, {x,y,z) € and {x,y,z) G M?. In the system 
of equations (1), a and jU can be seen as input parameters. 

To define the local system of coordinates, we have used a flat-Earth approach. 
As we want to analyse the motion of gUders during the TAEM phase, the height 
at which the TAEM phase starts is very small when compared to the Earth radius, 
justifying our analysis. However, this approach can be further refined by using an 
ellipsoidal coordinate system adequate to Earths shape, such as the WGS-84 coor- 
dinate system. 




Fig. 1 Local coordinate system for the point mass glider model. The origin of coordinates is lo- 
cated at the centre of mass of the aircraft, and the vector V is the velocity vector not necessarily 
colUnear with the aircraft longitudinal axis. 



In figure 2a)-b), we show the angle of attack a defined as the angle between 
the longitudinal reference line of the aircraft and the vector velocity of the aircraft. 
In airplanes, the angle of attack is always a positive angle. While in most aircrafts 
attack angles are always smaller than 15", the Space Shuttle is capable of attack 
angles up to 45", [8] and [7]. In figure 2c), we show the bank angle, defining the 
incUnation of the aircraft in the plane containing the velocity vector. 

The drag and lift forces in the system of equations (1) are given by, 

D{a,Ma) = qSCD{a,Ma) = \p{z)V^SCDia,Ma) 

(2) 

L{a,Ma) =qSCL{a,Ma) = ^p{z)V^SCL{a,Ma) 

where q = p {z)V^/2 is the dynamic pressure, S is the wing area of the aircraft, p (z) 
is the atmosphere density as a function of altitude (Appendix) and Ma is the Mach 
number. For each specific aircraft, the functions CD{(x,Ma) and CL{cc,Ma) are the 
aerodynamic drag and lift coefficients determined in wind tunnel experiments. 

Introducing the expressions (2) into equations (1), we obtain the final form for 
the equations of motion of a glider. 
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Fig. 2 In a) and b), we show the flight path angle 7 and the angle of attack a of an aircraft. In 
c) we show the bank angle \Jl, measuring the inclination of the aircraft in the plane containing the 
velocity vector and the horizontal direction. The flight path angle 7 depends on the angle of attack, 
on the aerodynamic coefficients of the aircraft and on the Mach number. The control of a glider is 
done by the manipulation of the angles of attack and bank. 



y = -^(z)sinr- iy^p{z)SCu{a,Ma)^ 



^ C0S7+ ^—p(z)5CL(a, Ma) ) VcoSjU ,< 
2m J cos 7 



' i = Vcosxcos 7 
j = y sin;ijcos7 (3) 
_ t = y sin7. 



where p (z) is calculated in the Appendix. 

The aircraft gliding trajectory is described by the system of equations (3), en- 
abling a simple geometric solution of the gliding aircraft control problem. 

When a glider is falling under a gravity field it converges to a steady state motion 
with a constant velocity and constant flight path angle given by, [1], 




1 



pS (C2+c2cos2;u)i/4 
Cp 
Clcosh 



(4) 



The geometry of the solutions of equation (3) in phase space are analysed in detail 
in [1]. 



3 Phenomenology of Space Shuttle gliding motion 

Using wind tunnel data for the operational range of aircrafts during the TAEM 
phase, we have done fits for the aerodynamic drag and lift coefficients Co and Cl of 
the Space Shuttle and these are well described by the parameterised functions, 
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CL{a,Ma) = (ai +a2a+a3a^)K{Ma)''i+"^2 
CD{a,Ma) = (0.01 +/iMa/2 +d3a^)K{May^+"'2 

where, 



K{Ma) = 




is a simplification of the Van Karman functions expanded to supersonic regimes, [8]. 
In table 1, we show, for the Space Shuttle, the parameter estimation of expressions 
(5) and (6) with wind tuimel data. 



Parameter 


Estimated 


Standard error 


t-statistics 


P-value 


fll 


-0.053 


0.009 


-6.15 


9.8 X 10"** 


02 


2.73 


0.06 


43.0 


1.8 X lO""-* 


03 


-1.55 


0.09 


-18.0 


2.0 X lO"^"* 


bi 


-1.01 


0.09 


-11.3 


7.4 X 10-"" 


bi 


1.1 


0.1 


8.7 


7.6 X 10-'^ 


ds 


1.79 


0.02 


99.0 


1.1 X lO-*"-^ 


ei 


-1.4 


0.1 


-12.6 


1.2x 10-" 


€2 


1.5 


0.1 


11.3 


5.8 X lO-"" 


/l 


0.028 


0.004 


6.46 


2.9 X IQ-'* 


/2 


1.4 


0.2 


8.57 


1.0 X 10-" 




1.25 


0.03 


49.6 


1.0 X 10-'"' 



Table 1 Parameters of the aerodynamic drag and lift coefficients (5) for the Space Shuttle, esti- 
mated from wind tunnel data, [6]. The significance of the fits have been determined with a chi- 

squared test. The large values of the absolute value of the t-statistics measures the likelihood of the 
parameters in the fits. The low values of the p-values mean that the fits are highly significant and 
the probability of finding a value outside the fitted ones are in the range IQ-* — IQ-^^. 



Introducing the expressions of CL{a,Ma) and CDiciC,Ma) into (4), changing the 
angle of attack a and the bank angle n leads to changes in the local steady states of 
the glider (see (4)), enabling a guided control of the direction of motion and of the 
glider speed. 

To control the aerodynamic behaviour of an aircraft, two main parameters are 
under the control of the aircraft commands: i) the bank angle n, and ii) the attack 
angle a. 

The bank angle n determines the inclination of the aircraft and is used for turn 
manoeuvres, figure 2c). 

The no-lift angle a„L, the max-glide angle a,„axgi and stall angle (Xstaii are partic- 
ular limits of the angle of attack of an aircraft, figure 2a) and 3. 

In figure 3, we show, for several values of the Mach number, the behaviour of the 
ratio L/D, as a function of the angle of attack a for the Space Shuttle. All the curves 
intersect at the no-lift angle a„L. The no-lift angle a„L is the angle for which L/D 
is zero due the absence of the lift force and is independent of the speed. The max- 
ghde angle OCmaxgi is the angle that maximises the ratio L/D, and is dependent on 
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the Mach number. The stall angle CCgtaii is the angle at which hft dependency looses 
linearity and lift peaks before beginning to decrease. The stall angle is independent 
of the Mach number. 

6 
5 
4 

d 2 
1 


10 20 30 40 50 
a (degrees) 

Fig. 3 Ratio L/D, as a function of the angle of attach a, for the Space Shuttle at different Mach 
numbers, calculated from (5)-(6) and table 1. The no-lift parameter is a„L = 1.5", the stall angle is 
CCstall = 45° and a,„axgl is given by (7). The Space Shuttle is a glider and thus can only move across 
its L/D curve. For higher speeds this curve will become increasingly flat and the max-glide angle 
Omaxgl will move further to the right reaching saturation. 

With the functions (5)-(6), we have approximated the max-glide angle OCmaxgi as 
a fimction of Mach number. For the case of the Space Shuttle, we have obtained, 

_ f 0.0906 + 0.0573Ma + 0.007 lMa2 {Ma < 1.25) 

~ \ 0.1070 + 0.0577Ma-0.0037Ma2 (1.25<Ma<5) *■ 

determined with a correlation coefficient of = 0.999. The Mach number is defined 
by Ma = Y sound where the sound speed is calculated with, 

%ound = VWWs (8) 

and r(z) is given in table 2 in the Appendix. 7= 1.4 is the diatomic gas constant 
and/?, = 287.04 J/(kgK). 

4 Dynamic trajectory control of gliders 

A gUder is not always in an equilibrium state but naturally converges to it given 
enough time. Our algorithm will take advantage of this behaviour by determining 
the equilibrium conditions needed to reach the target, imposing them on the system 
and letting it evolve in time. 

To define the control problem, we consider the initial condition. 



Ma=0.6 








f ffstall 
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(xo,yo,2o,Vo,ro,Xo) 
defining tlie initial coordinates of the TAEM phase. Let, 

be the space coordinates of the target, which coincide with the central point in the 
HAC region. We consider that the target point is only defined by the spatial coordi- 
nates of the HAC, and the direction of the velocity vector is arbitrary. In fact, this 
is possible at low altitudes (3 km) because the atmosphere is dense enough to allow 
the ghder to preform turns in short distances and the vehicle is always travelling 
near the equilibrium speed. 

The intermediate coordinates of the glider path are, 

i^i jyijZi,Vi,Yi, Xi ) 

where i = 0, These intermediate coordinates are evaluated at time intervals 
T 

con ■ 

In the configuration space {x,y, z), we define the direction vector from the current 
position of the glider to the target point as, 

Pi = {^f - , yf - yi ,zf-zi). (9) 

In order to direct the aircraft to the target, we control the attack and bank angles 
separately. 

In the attack angle heading control, we analyse the glider trajectory in the three 
dimensional ambient space (x,>', z), and we command the gUder trajectory path angle 
by controlling the angle of attack a. 

In the bank angle heading control, the control procedures will be done in the 

{x,y) plane by adjusting the bank angle jU. 

At the step number / of the dynamic control process, the initial conditions are 
{xi,yi,Zi,Vi,Yi,Xi)- At this stage, the angle of attack and bank angle are a, and /i,-. 
Then, we calculate the new values of the glider control parameters a,+ i and i by 
the two procedures described below. With these new values for a and jU , the aircraft 
will follow a new trajectory during the time interval Tcm, figure 4. 

This control process is done sequentially in time, until the glider reaches de HAC 
region. In practical terms, the control mechanisms stops when the distance from the 
spacecraft to the centre of the HAC point attains a minimum. 

We analyse now in detail the two control and conomand procedures for a and /i. 

Attack angle heading control: 

The attack angle heading command and control was designed so that the vehicle 
is always re-orienting vertically to the HAC point through a straight Une path. 
The tangent of the angle between the x — y projection and the z component of 
the direction vector Pi to the target point is computed at each iteration, and we 
obtain. 
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Start 


Measurements: 

Time (t), Position, Velocity 

vector and atmosplieric state 

variables. 










Calculate new attack angle and new 
bank angle. 

Change the glider controls to the new 
values. 




i 




Wait seconds 


No 


<Cr HAC target height reachedj>- 



Yes 






END 





Fig. 4 Block diagram for the controller algorithm. 



Gi+i 



^{xf-Xiy + iyf-yiY 



where {xi,yi,Zi) is the current position of the gUder. At this position, the gUder 
has flight path ji- Then, to direct the motion of the ghder to the target with a 
steady flight path, by (4), we must have. 



Gi+i =tan7 = - 



1 



Cl /Cd cos Hi+i 



(10) 



Assuming that it is possible to direct the motion to the target using a nuU bank 
angle, ;U,+i = 0, we solve equation (10) in order to the ratio Cl/Cd, and we obtain 
the solution c/+i. Then: 

a) If c,+i is bigger flian CL^OmaxghMai) /CoicimaxghMai), the target cannot 
be reached in a straight-line and the max-glide attack angle will be selected, 
ai+i = oCmaxgi- The curve of Cl/Cd as a function of a and of the Mach number 
Ma is given by (5) and (6), and (Xmaxgi is calculated from (7) and (8). 

b) If a+i is smaller than CL{oCstaii,Mai)/CD{oCstaii,Mai), the target cannot be 
reached in a straight-Une and the stall angle will be selected, a,+i = oCstaii- 

c) Otherwise, the attack angle ttj+i is computed by solving the equation 
CL{a,Mai)/CD{a,Mai) = Cf+i. 

At this stage, we have chosen a new attack angle o;,+i . With this new attack angle, 
we re-orient dynamically and vertically the aircraft trajectory to the target. 



Bank angle heading control: 
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The bank angle heading control was constructed in such a way that, in the {x,y) 
plan, the aircraft is always re-orienting horizontally to the HAC. 
The angular misalignment between the direction vector to the target point (9) and 
the speed in the {x,y) plane is measured using the dot product. The direction is 
measured by the j component of the exterior product (A) between the direction 
vector to the target point P, and the aircraft speed Vj. With P- = Pi^ex + PiySy and 
V; = Vi^ex + Vi^ey, in order to align the aircraft to the target point in the (x, >•) 
plane, the new bank angle is. 




P' v 

' Sign((P, AV,).) 



P<l|x||V; 



(11) 



where, we have introduced a new constant Thard € [0) 1] ■ The higher this constant, 

the faster the vehicle will turn for the same angular deviation. 

We impose now a security threshold in the bank angle, Umax- A typical value for 

the maximum bank angle is jj.max = ±70". Therefore, the new control bank angle 

is, 

= min{|Mf_^f|, |Mm«.|}.Sign(Mf+1) • (12) 



5 Simulations 

In the previous section, we have described a control mechanism in order to guide a 
gUder to a target. At each time step, the algorithm determines the shortest path to 
the target and determines the unique values of the attitude commands of the glider 
that are compatible with the aerodynamic characteristics of the glider. We now test 
this algorithm with some numerical simulations. 

We have taken the glider initial coordinates (xo,>'o,zo) = (0,0,40000) m, Vo = 
1000 m/s, Yo = 0, xo=0, ^0 = 0", ao = 30", T^on = 0.1 and Thard = 1-0, and we 
calculated the trajectories of the glider by numerically intreating equations (3) with 
a fourth order Runge-Kutta integration method. 

The goal was to reach some target point that we have defined as the centre point 
of the HAC. We have chosen three different target HAC points with coordinates, 

1) (x/,y/,z/) = (200000, 10000,3000) m (figure 5). 

2) (x/,>'/,z/) = (50000, 10000,3000) m (figure 6). 

3) {xf,yf,zf) = (0, 10000,3000) m (figure 7). 

and we have calculated the controlled trajectories from the same initial point. The 
arrival to the HAC point occurs when the distance from the glider to the centre of 
the HAC point attains a minimum. This distance error will be denoted by e,/. In 
figures 5, 6 and 7, we show the glider controlled trajectories as function of time and 
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Fig. 5 Trajectory of the glider in the ambient space and control commands as a function of time. 
The coordinates of the HAC point target are {xf,yf,Zf) = (200000, 10000,3000) m. The time of 
arrival at the HAC isf = 539.6 s, with a distance error ej = 14.6 m and final speed V/ = 0.203 M. 
The dots indicate the position of the glider after 100 s and 200 s of flight and the HAC position. 
Positive values of correspond to left turns and negative values of /J correspond to right turns. 



the sequence of the attack and bank angle values as computed by the command and 
control algorithm. We have computed the time of arrival at the HAC, the final speed 
at the HAC (Vf) measured in Mach number units, and the distance error e^. 

The basic features of this algorithm is to guide the aircraft to the HAC point 
with very low distance errors. The choice of the initial conditions has been done 
insuring that the initial energy of the ghder is enough to arrive at the target point. 
In this study, we have chosen target points within the maximum range calculated 
numerically by imposing the condition that the flight is always done with zero bank 
angle and maximum ghde angle. In this case, the ratio Cl/Cd is maximal and the 
drag on the glider is minimal. For the initial conditions chosen and the Space Shuttle 
parameters, the range is of the order of 286 km. 

Dynamic aircraft trajectories computed with the algorithm presented here depend 
on the control time Tcon- For the conditions in figure 5, we have evaluated the dis- 
tance error from the centre of the HAC as a function of Tcon- For ^con < 30, we have 
found that, 

= 13.7e'' '^9^- . (13) 

In figure 8, we show the dependence of the distance error on the control time T^on 
for the initial and final conditions of the simulation in figure 5. 
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Fig. 6 Trajectory of the glider in the ambient space and control commands as a function of time. 
The coordmates of the HAC point target are (x/,y/,z/) = (50000, 10000,3000) m. The time of 
arrival at the HAC isf = 345.9 s, with a distance error = 23.1 m and final speed Vf = 0.205 M. 
The dots indicate the position of the glider after 100 s and 200 s of flight and the HAC position. 




Fig. 7 Trajectory of the glider in the ambient space and control commands as a function of time. 
The coordinates of the HAC point target are (xf,yf,zf) = (0, 10000,3000) m. The time of arrival 
at the HAC t = 485.9 s, with a distance error e,; = 51.3 m and final speed Vf = 0.200 M. The dots 
indicate the position of the gUder after 100 s and 200 s of flight and the HAC position. 
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Fig. 8 Dependence of the distance error on the control time Tcon for the initial and final conditions 
of the simulation in figure 5. For Tcon < 30, the distance arrow follows the approximate exponential 
law (13), represented by the thin Une. 



We have also tested the dependence of the controlled trajectories as a function of 
the entry angle Xq. In figure 9, we show the trajectories as in figure 5 but with Xa = 
;r/4,0, —n/A. In this three cases, the distance errors are = 34.6 m, = 14.6 m 
and Cd = 52.2 m, respectively. For larger values of angles Xo, the distance error can 
be as large as 69 km (Xo = 7t/2). 




50000 100000 150000 200000 
x(m) 



Fig. 9 Different trajectories calculated with the initial parameter Xo = Jr/4,0, — ;r/4. The other 
parameters are the same as in figure 5. For these trajectories, the distance errors are ej = 34.6 m, 
ed = 14.6 m and = 52.2 m, respectively. 



6 Conclusions 

We have derived a new algorithm for the command and control during the TAEM 
phase of re-usable space vehicles. The algorithm determines locally the shortest path 
to the target point, compatible with the aerodynamic characteristics of the aircraft. 
We have tested the ability of the algorithm to guide the Space Shuttle during the 
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TAEM re-entry orbit, proving the feasibility of the algorithm, even using control 
times of the order of 30 s. Further refinements of the algorithm are under study [1]. 
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Appendix 

The Earth atmosphere parameters are based on the 1976 US Standard Atmosphere 
Model. For the first seven layers we have used the formulas described in [10]. In 
table 2, we show the parameterisation of the thermodynamic quantities for the Earth 
atmosphere. 



Layer 


zo (m) 


To (K) 


A(,(K/m) 


Po (Pa) 


1 





288.15 


-0.0065 


101325.00 


2 


11019 


216.65 




22632.10 


3 


20063 


216.65 


0.0010 


5474.89 


4 


32162 


228.65 


0.0028 


868.02 


5 


47359 


270.65 




110.91 


6 


51412 


270.65 


0.0028 


66.94 


7 


71802 


214.65 


-0.0020 


3.96 



Layer 


r(K) 


P(Pa) 


P (kg/m^) 


1 


7b + Ao(z-zo) 


Pg(5l)«(z)M<„v/{ifA„) 


p 

TR, 


2 


To 




P 
TR, 


3 


7b + Ao(z-zo) 




P 

TR, 


4 


To-\-^{z-zo) 


Pg(a)s(z)M<„v/(;?Ao) 


P 
TR, 


5 


To 




P 
TR, 


6 


7b + Ao(z-zo) 




P 
TR, 


7 


To + ^{z-zo) 


PQ('A)g{malr/{Rh) 


P 

TR, 



Table 2 Characteristic parameters for the lower layers of the atmosphere, zo is the lower alti- 
tude of the layer, R = 8.31432 J/(mol kg) and R, = 287.04 J/(kg K) are gas constants, Ma,> = 
0.0289644 kg/mol, g(z) = 8o{Re / {Re + z)f is the gravity acceleration, go = 9.80665 m/s^ is the 
standard gravitational acceleration constant and Re = 6.371 x 10* m is the Earth mean radius. 
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